Reduced density matrix functional theory at finite temperature. I. 

Theoretical foundations 
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An ab-initio approach to the description of grand canonical ensembles in thermal equilibrium, 
having local or nonlocal external potentials, will be presented. To this end a variational principle 
for the grand potential ft of the system under consideration with respect to its one-reduced density 
matrix 7 will be established. The domain of ^[7] will be shown to be determined by a simple set of 
constraints, making it suitable for numerical minimization. We will furthermore prove the existence 
of an analogon to the Kohn-Sham system of density functional theory, i.e. a system of noninteracting 
particles reproducing the one-reduced density matrix of the interacting system at finite temperature. 
Starting from this Kohn-Sham system as unperturbed system, we deduce a many-body approach to 
iteratively construct approximate functionals for the grand potential. 



I. INTRODUCTION 

Density functional theory (DFT), based on the cele- 
brated theorems of Hohenberg and Kohn [1], is today's 
method of choice for the description of groundstate (gs) 
properties of quantum systems. There are, however, 
cases in which DFT performs rather poorly. Prominent 
examples are the so-called Mott insulators, for which 
some of the transition metal oxides constitute the best- 
studied representatives. For this type of strongly cor- 
related systems, the Kohn-Sham description with stan- 
dard exchange-correlation functionals typically yields a 
metallic ground state while the true experimental low- 
temperature phase is that of an antiferromagnetic insu- 
lator. 

At finite temperature the description of strongly cor- 
related systems is even more challenging. Genuine Mott 
insulators are characterized by the feature that, when 
heated up from their antiferromagnetic insulating gs, 
they stay insulating above the Neel temperature, i.e. in 
the absence of long-range magnetic order. By contrast, 
weakly correlated "Slater" insulators become metallic at 
the Neel temperature. 

One possible way to tackle this challenge is to search 
for more accurate functionals in the framework of DFT. 
Alternatively one may look for other theoretical frame- 
works in which the treatment of strong correlation might 
be easier. One candidate for such a framework is reduced 
density matrix functional theory (RDMFT). Through its 
more direct treatment of many-particle correlations it has 
shown big promise for the description of finite [2-7] as 
well as infinite systems [8-10]. In particular, it was pos- 
sible to predict insulating ground states for transition 
metal oxides without breaking the spin symmetry [10]. 

Motivated by the success at zero temperature, the pur- 
pose of the present work is to lay the theoretical founda- 
tion for the finite-temperature version of RDMFT (FT- 
RDMFT). As this theory will be a truly general ab-initio 
theory its applicability is not restricted to the case of 
Mott insulators. There is a wide variety of physical phe- 



nomena whose investigation require an accurate descrip- 
tion of temperature effects. These phenomena include 
temperature driven magnetic [11, 12] or superconducting 
[13, 14] phase transitions in solids, properties of shock 
compressed noble gases [15, 16], femtochemistry at sur- 
faces of solids [17], the properties of plasmas [18-20], and 
the composition of different phases of water in the inte- 
rior of giant planets [21]. 



II. OUTLINE 

To lay the foundation for FT-RDMFT we will show 
that at finite temperatures the grand potential of sys- 
tems with generally nonlocal external potentials can be 
written as a functional of the one-reduced density matrix 
(1RDM). The domain of the functional is then shown to 
be determined by just a few simple constraints. 

We will prove the existence of a one-to-one mapping 
between the 1RDM at equilibrium and the correspond- 
ing statistical density operator (SDO). Furthermore, we 
will show that at finite temperatures there exists a nonin- 
teracting Kohn-Sham system with possibly nonlocal ex- 
ternal potential which reproduces the 1RDM of the fully 
interacting many-body system. 

The above mentioned facts have several important con- 
sequences. Firstly, the Kohn-Sham system can be em- 
ployed to establish self-consistent minimization proce- 
dures which is not possible at zero temperature (where 
no Kohn-Sham-like system exists) [22, 23]. Secondly, 
the existence of this Kohn-Sham-like system also has the 
benefit of giving easy access to one-electron properties 
(like in DFT). Thirdly, the domain of external poten- 
tials includes nonlocal potentials (e.g. pscudopotcntials 
[24, 25]), for which the original Hohcnbcrg-Kohn theo- 
rems of DFT and FT-DFT do not hold [26]. 

Furthermore, based on the Kohn-Sham system of FT- 
RDMFT we shall derive a methodology to iteratively con- 
struct functionals of the 1RDM by employing finite tem- 
perature many body perturbation theory (FT-MBPT). 
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We will conclude this work by considering the case of 
fixed particle number TV, i.e. we will discuss how the re- 
sults for grand canonical ensembles differ from the ones 
for canonical ensembles. 

Part II [27] and III [28] of this work will investigate the 
performance of FT-RDMFT when applied to the homo- 
geneous electron gas. Part II will focuss on the first- 
order, i.e. exchange-only, functional considering both 
collinear as well as spin spiral states. Part III will then 
discuss correlation contributions in the framework of FT- 
RDMFT focussing on correlation energies, momentum 
distributions and phase diagrams for collinear spin con- 
figurations. 



III. FINITE-TEMPERATURE RDMFT 

We start from the standard definition of the SDO D, 
which is a weighted sum of projection operators on the 
underlying Hilbcrt space. 



> 0, 



i. 



(i) 



If one considers grand canonical ensembles, one has to 
allow for a change of particle number. The appropriate 
Hilbert space is then given as the direct sum of sym- 
metrized tensor products of the one-particle Hilbcrt space 
h which is called the Fock space W, 



H = Q)Sh 



(gm 



(2) 



As long as the system does not allow a mixing of states 
with different particle numbers it is possible to reduce the 
set of all possible SDO from Eq. (1) and only use pro- 
jection operators on states with defined particle number 
N. 

b 9C = ^2w aN \y aN )(y aN \, 

a,N 

WaN >0,/] W aN = 1 (3) 

The main thermodynamic variable in a grand canonical 
ensemble is the grand potential il. For a given SDO ft is 
defined as 



n[t>] = tr{D(H - fxN + 1/pinD)}, 



(4) 



where H is the Hamiltonian of the system. The two 
Lagrangian multipliers /i (the chemical potential) and 
1//3 (the temperature) describe the coupling to the par- 
ticle and heat baths. The thermal equilibrium of a grand 
canonical ensemble is then defined as that SDO for which 
the grand potential Q,[D] is minimal. This definition of 
the equilibrium of a grand canonical ensemble leads to 



the finite-temperature Rayleigh-Ritz variational princi- 
ple (the Gibbs principle) [29] 



with 



si[D]>n[b eq \, b^b eq , (5) 



D eq = e-^-^/tT{e-^-^}. (6) 



The 1RDM 7(2;, x'), corresponding to the SDO D, is de- 
fined as 



70,a/) = tr{DiP+(x')^{x)} 



(7) 



The variable x denotes a combination of spin index a and 
spacial coordinate r (x = (r, a)). An integration over x is 
therefore to be interpreted as an integration over r and a 
summation over a. By construction, 7(2:, x') is hermitean 
and can therefore be written in spectral representation as 



(8) 



with real- valued eigenvalues {^i}, which are called occu- 
pation numbers, and eigenfunctions {4>i(x}}, which are 
called natural orbitals, [30]. The neccessary and suffi- 
cient conditions for N-representability of j(x, x') are [31] 
that {4>i} is a complete set and 



< Hi < 1, ^rii = N 



(9) 



It is sometimes desirable to treat spin and spacial vari- 
ables seperately. To this end we introduce a two- 
component (Pauli) spinor notation. 



Mr) 



0ii(r) 
<Mr) 



(10) 



where <pi a (r) = <fii(x) = <pi(a,r) (a = 1,2) are the or- 
bitals of Eq. (8). The 1RDM can then be written as a 
matrix in spin space as 
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(r,r') =^ni$j(r')(g)$i(r) 



0Mr%i(r) ^(r')^i(r 
#!(r% 2 (r) ^(r')fe(r 



(11) 



(12) 



There are some cases in which one can treat different 
spin channels seperately. One of these special cases is 
the particulary important situation of collinear spin con- 
figuration. For these systems, the natural orbitals are 
so-called spin orbitals, i.e. spinors containing only one 
spin component. We characterize these spin orbitals by 
an additional index. 



*<i(r) 







,$ i2 (r) 





<fe(r) 



(13) 



This then leads to a 1RDM where every 2x2 matrix in 
Eq. (12) contains only one nonvanishing entry, either the 
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11 or the 22 one. Hence the complete 1RDM is diagonal 
w.r.t. the spin coordinate, 

7 CTCT '(r,r') = <WX^C( r ')<Mr), (14) 

i 

where Ti^fj &.re the occupation numbers of the special 
spinors <&i a {v) of Eq. (13). Another special case where 
this seperation is possible is the case of spin spiral states 
which will be dealt with in Part II [27] of this work. 

We will now turn to the question of the existence of 
a Kohn-Sham system in RDMFT, i.e. a nonintcracting 
system which reproduces the gs-lRDM or eq-lRDM of 
an interacting system. When dealing with Coulomb sys- 
tems, i.e. systems where the two-particle interaction is 
w(x,x') = l/|r — r'l, the following theorem was shown 
[32]. 

The gs-lRDM of a Coulomb system at zero temper- 
ature has an infinite number of occupied orbitals, i.e. 
natural orbitals with occupation numbers > 0. 

The original proof is rather involved and we therefore 
present an alternative simpler version in the following. 

As Kato showed in 1957 [33], the divergence of the 
Coulomb interaction for vanishing inter-particle dis- 
tances leads to a nucleus-electron as well as an electron- 
electron cusp in the eigenfunctions of the Hamiltonian. 
Every eigenfunction is an element of the general N- 
particle Hilbert space of the system and can therefore 
be expanded in any basis of this Hilbert space. Further- 
more, we know that one can create an N-particle ba- 
sis by constructing all possible N-particle Slater determi- 
nants of a given one-particle basis. As in a Slater deter- 
minant a particular one-particle wavefunction "does not 
know" about the coordinates of another one, the electron- 
electron cusp condition cannot be fulfilled by a super- 
position of a finite number of Slater determinants. One 
therefore needs an infinite number of Slater determinants 
to reproduce an eigenfunction of a Coulomb system. This 
argument is valid for all one-particle bases, including the 
set of natural orbitals of the gs-lRDM. If an infinite 
number of N-particle wavefunctions contributes, then the 
number of contributing one-particle natural orbitals also 
has to be infinite, implying the statement above. 

As the ground state wavefunction of a nonintcracting 
nondegenerate system is a single Slater determinant, the 
gs-lRDM of a nonintcracting system will be idempotent. 
Therefore, we conclude that at zero temperature there 
exists no Kohn-Sham system for Coulomb systems. We 
would like to mention, however, that a statistical mix- 
ture of the degenerate ground states of a nonintcracting 
system is capable of reproducing a 1RDM with partially 
occupied natural orbitals [34] but the degree of degener- 
acy is equal to the number of partially occupied states. 
For a Coulomb system therefore the degeneracy must be 
infinite. 

We will show in the following that at finite temperature 
this disadvantage disappears, by proving the existence 
of a Kohn-Sham system in FT-RDMFT. We will show 



how the Kohn-Sham Hamiltonian can be constructed di- 
rectly, i.e. from an explicit analytical formula, from the 
occupation numbers and natural orbitals of the interact- 
ing 1RDM. This direct construction is an advantage over 
DFT where, if one wants to evaluate the Kohn-Sham po- 
tential for a given density, the single-particle Schrddinger 
equation must be inverted numerically. This is in gen- 
eral a nontrivial task which can be achieved by iterative 
methods [35]. 



A. Kohn-Sham system for FT-RDMFT 

In order to show the existence of a Kohn-Sham system 
in the context of FT-RDMFT, it is instructive to consider 
an arbitrary noninteracting fcrmionic system defined by 
the one-particle Hamiltonian Hs with eigenvalues {e^} 
and eigenfunctions {(f>i(x)}. 

6™ =£>|&)<&| (15) 

i 

For a grand canonical ensemble, the eq-lRDM is then 
[36] given by 

-f eq (x,x') = ^2m<l)*(x')(j)i(x), (16) 

i 

where the occupation numbers {n,} are determined com- 
pletely by the eigenvalues {ei} and the chemical potential 

This relation can be inverted to give 

'■-"'H^r)- (18) 

From Eq.(18), it is now possible to construct the Kohn- 
Sham Hamiltonian from a given 1RDM. Its eigenfunc- 
tions are given by the natural orbitals of the 1RDM while 
the eigenvalues are defined by Eq.(18). As can be seen 
from Eq. (17) the occupation numbers of the eq-lRDM 
of a nonintcracting system in thermal equilibrium at fi- 
nite temperature cannot be or 1. To ensure, that the 
domain of the noninteracting grand potential functional 
will include the domain of the interacting one we show 
in the Appendix A that the grand canonical eq-lRDM of 
an interacting system will also have occupation numbers 
strictly between and 1. 

Defining the kinetic operator in the basis of natural or- 
bitals tij = ((pi\T\4>j), the effective onc-particlc potential 
v e ff{x,x') can be expressed explicitly as 

v eff {x,x') = ^J<V , - Uj^ix'^jix). (19) 

i,3 
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As we can see from Eq. (19), the noninteracting poten- 
tial which reproduces a given 1RDM is generally nonlo- 
cal. This shall serve as motivation to prove a Hohcnbcrg- 
Kohn theorem for systems in grand canonical equilibrium 
subject to arbitrary interactions and nonlocal external 
one-particle potentials. 

B. Hohenberg-Kohn theorems for FT-RDMFT 

For local (multiplicative) external potentials the 
Hohenberg-Kohn theorem (at zero temperature) and the 
Mcrmin theorem (at finite temperature) immediately im- 
ply a one-to-one correspondence between the potential 
and the 1RDM. This can be seen from the following re- 
lations. 

7(r,r') — >p(r) v(r) (20a) 
v(r) — ► *(£>) — ► 7 (r, r') (20b) 

Every 1RDM yields the density as its diagonal and the 
Hohenberg-Kohn or Mermin theorems prove the one-to- 
one correspondence between the density and the poten- 
tial (see Eq. (20a)). This potential, on the other hand, 
determines the Hamiltonian and therefore the ground- 
state \Jj as well as the eq-SDO D which in turn lead to 
the 1RDM via Eq. (7) (at zero temperature for a non- 
degenerate system the SDO is given as a single projector 
on the groundstate) . Eqs. (20a) and (20b) are therefore 
closed and the Hohenberg-Kohn and Mermin theorems 
imply one-to-one correspondences between the 1RDM, 
the density, the potential and the groundstate, i.e. the 
eq-SDO. 

In 1974 Gilbert [26] extended the zero-temperature 
Hohenberg-Kohn theorems to systems with nonlocal ex- 
ternal potential. He showed that for such systems the 
ground state is not uniquely determined by the density 
anymore and one has to resort to the 1RDM. 

In the following, we will show that at finite tempera- 
ture the 1RDM is still sufficient to describe the equilib- 
rium properties of an interacting system with nonlocal 
external potential. 

The three ingredients needed to establish the founda- 
tion for FT-RDMFT are: (I) the map between D eq and 
r y eq (x,x') has to be invertible, implying the existence of 
a grand potential functional Sl[y], (II) the existence of 
a universal internal-energy functional F[y] and (III) the 
minimization of ^[7] leads to the eq-lRDM. 

(I) The proof of the existence of a one-to-one map- 
ping between the eq-SDO D eq (see Eq.(6)) and the eq- 
1RDM (see Eq.(7)) can be divided in two parts. First, 
the one-to-one mapping between D eq and the external 
potential minus the chemical potential (v ex t(x,x r ) — fx) 
can be shown, then the one-to-one correspondence be- 
tween (v ex t(x,x') — /x) and j eq (x,x') is proven. 

•D eq «-» (v ext (x,x') - fx) 
Let H and H' be two different Hamiltonians and assume 
they lead to the same SDO D. H' shall differ from H 



only by a one-particle potential contribution U . With 
Eq. (6) this reads 

e - Pi H-^N) /z = e -p { H+U-»N) /z ,_ (21) 

where Z and Z 1 are the partition functions (e.g. Z = 
tj .{ e -p(H-»N)yy Solving Eq. (21) for U yields 

U = f dxdx'u{x' ,x)^ + {x')^{x) = -~ln— , (22) 
J P z 

We now need to show that there is no one-particle po- 
tential u(x, x') 7^ which fulfills this equality, thereby 
contradicting our initial assumption. To proceed, we 
assume three different Slater determinants \Xi) = 
|1,0,0,...>, |X 2 ) = |0,1,0,...>, and \X 3 ) = |1, 1,0, ...> in 
the basis {xi}- The potential in this basis is denoted 
by Ui 3 ; = J dxdx'u{x', x)x*(x')xj(x). Calculating the ex- 
pectation value of both sides of Eq. (22) w.r.t. these 
three Slater determinants we get the following system of 
equations 

1 Z 

-3 In— = uu = U22 = u n +u 2 2, (23) 
P 6 

which can only be fulfilled by u\\ = U22 = and Z = Z' . 
A repetition of this argument for all possible bases then 
shows that only U = fulfills Eq. (22) which in turn 
proves the one-to-one correspondence between D eq and 
(v ex t{x, x') —fi). It has to be noted that this proof intrin- 
sically relies on the fact that in the case of grand canon- 
ical ensembles we work in a Fock space, i.e. a Hilbert 
space with varying particle number. If we had restricted 
ourselves to a canonical situation, i.e. a Hilbert space 
with fixed particle number, we would have found that 
the external potential is uniquely determined only up to 
an additional constant. 

This proof is valid for any finite temperature. It is 
based on the bijectivity of the exponential function which 
allows us to invert Eq. (21), leading to Eq. (23). At zero 
temperature this bijectivity breaks down and, in fact, 
a general one-to-one correspondence (v ex t (x, x') — n) 
\I/ cannot exist which is easily shown by the following 
counterexample. Assume the following two potentials 

v(x, x') = vhf{x, x') (24) 
v'(x, x') = v H f{x, x') + a~/ gs {x, x'), (25) 

where Vhf{x, x') denotes the Hartree-Fock potential and 
7 gs (a;, x') is the corresponding gs-lRDM. Both potentials 
yield the same gs-dcterminant for systems with a discrete 
spectrum and sufficiently small a. We were able to show 
in Appendix B that in general, the one-to-one correspon- 
dence between potential and wavefunction persists if and 
only if there are no pinned occupation numbers, i.e. no 
occupation numbers equal to or 1. The case described 
above exhibits only pinned states and can therefore be 
be seen as an extreme case of our general result from 
Appendix B. 
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• (v ext (x,x') - fi) o i{x,x') 
In order to prove the one-to-one correspondence between 
(v ex t{x, x') — [i) and 7 eg (a;, x'), consider H and ii 7 differ- 
ing only in their external and chemical potentials. The 
corresponding grand potentials are given by 

Q[D eq ] = tr{D eq (H -{iN+ l/p\nD eq )} (26) 
n'[D' eq \ = tr{I)' eq (H' - n'N + l//31n^,)}, (27) 

where D eq and D' eq are defined according to Eq. (6). As 



from the first part of our proof we know that D eq ^ D' eq 
the variational principle (Eq. (5)) leads to 

Q[b eq ] < n[b' eq ] (28) 

= tr{D' eq (H-ixN + l/p]nD' eq )} (29) 

+ tr{& eg ((H-ixN)-(H' -n'N))} (30) 
By exchanging primed and unprimed objects one obtains 



Q[D eq ] < n'[D' eq ] + J dxdx'{{v ext {x',x)-n) - (v' ext (x',x)~ f i')h'(x,x / ) 
n'[b' eq \<n[b eq ] + J dxdx'((v' EXt (x',x)- ti')-(v ext (x' : x)~ [i)h(x,x') 
Adding these two equations leads to the following relation; 

dxdx'((v' ext (x',x) -//) - (v e xt(x',x) - n))(i(x, x') -f'(x,x')) > 0. 

I 



(31) 
(32) 

(33) 



The existence of two different sets of external and chem- 
ical potentials yielding the same eq-lRDM lets the in- 
tegral in Eq.(33) vanish which leads to a contradiction. 
Hence the initial assumption is falsified. 

This proof of the existence of a one-to-one mapping 
between D eq and j eq is necessary and sufficient to define 
the grand potential as a functional of the 1RDM. 

n[ leq ] =tr{D[ leq ](H - »N +1/J3]n(b[ leq }))} (34) 

(II) In Eq. (34) the contributions from the exter- 
nal and the chemical potential can be separated yielding 
the definition of a universal functional -^"[7 ei j] for FT- 
RDMFT. 



^lleq] = J dxdx'{v e 



(x, x') - fji)j e g(x, x') + 7"[7, 



T[ leq ] = tr{D[ leq } (T + W + 1//3 In D[ leq }) } 



(35) 



The functional J-[y] is defined on the set of lRDMs com- 
ing from SDOs of the form of Eq. (6). The conditions 
to ensure that an arbitrary 1RDM comes from such a 
SDO arc unknown. This set of lRDMs is a subset of the 
set of all ensemble-N-representable lRDMs. Therefore, 
following ideas of Valone [37] and Lieb [38] , the domain 
of f2[7] can be extended to the whole set of ensemble- 
N-representable lRDMs. The new universal functional 
3^7] as well as the grand potential functional £l\y\ then 
become 

^[7] = J dxdx'(v ext {x,x') — fi)i(x,x') +#[7] (36) 
£[7] = inf tT{D(f + W+l/phib)}. (37) 



(The corresonding zero-tempcraturc formulation is also 
sometimes referred to as the grand canonical ensemble 
formulation [39]). 

(Ill) The variational principle in Eq. (5) now allows 
the determination of the equilibrium grand potential by 
minimizing the functional ^[7]. 



< 9 



min SIM 

7 er« 



(38) 



T N is the set of all ensemble-N-representable lRDMs. We 
would now like to postulate the following FT-RDMFT 
version of the Eulcr-Lagrangc equation for the eq-lRDM 

leq 



5j(x, x') 



+ v ext (x,x') = /i 



(39) 



7=9 



but there are some sublcties to be mentioned. Firstly, 
the question of existence of a functional derivative of the 
exact ^[7] at the eq-lRDM has not been answered so 
far and is subject to continued studies (also in DFT this 
question is still under discussion [40]). In an approxi- 
mate treatment, however, one will most likely define a 
functional F a PP rox [^/] for which the functional derivative 
exists, justifying the use of the Euler-Lagrangc equation. 
Secondly, if the eq-lRDM will be on the boundary of 
the domain T N , the Euler-Lagrangc equation does not 
take the form of Eq. (39) but has to incorporate the 
constraints on the eigenvalues of the 1RDM by further 
Kuhn- Tucker multipliers [41]. Fortunately, as we have 
shown in the appendix A, in the case of grand canonical 
equilibriums the minimizing 1RDM will never be on the 
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boundary of T N , rendering these additional multipliers 
unnecessary. 

We have shown that FT-RDMFT is theoretically well 
founded and that it can be used to describe general 
quantum systems in thermal grand canonical equilib- 
rium. This has a wide variety of possible applications 
[15, 18, 42-44] and it raises the hope that FT-RDMFT 
performs better than FT-DFT which has been success- 
fully applied to a rather limited class of problems [18, 19]. 

The central problem in FT-RDMFT now consists in 
the construction of correlation functionals. In the follow- 
ing we will exploit the existence of a Kohn-Sham system 
to derive a methodology for the iterative construction of 
functionals using methods from FT-MBPT. 



Construction of Functionals 



Following the construction of the adiabatic connection in 
DFT, we introduce a coupling constant A and define the 
following function 



0(A) 
Tr 



r{£>(A) (f + XW + V(\) + l/p\nD(X)^Y (46) 



where V(X) is chosen such that the eq-lRDM does not 
change on changing A. For full coupling strength it re- 
produces the physical external and chemical potentials 
whereas for zero coupling it represents the Kohn-Sham 
contributions. 



V(l) = V ext - [iN 
V(0) = V KS - msN 



(47) 
(48) 



Following ideas of DFT, we would like to express the 
functional for the grand potential of an interacting sys- 
tem as the grand potential of a noninteracting system 
with additional terms originating from the interaction. 

o[ 7 ] = n fc [ 7 ] + v ext [ 7 ] - m tv[ 7 ] - i/£So[7]+ 

Slwb] (40) 

The functionals of kinetic energy £^[7], external poten- 
tial T4xt[7], particle number iV[ 7 ], and noninteracting en- 
tropy So [7] are explicitly known. 



fife [7] = 


1 dx' lim f - ^-) l( x > 


X>) 


(41) 


v ext M = 


j dxdx'v ex t(x' , x)j(x, x 


') 


(42) 


AT[ 7 ] = 


j dx^f(x,x) 




(43) 


So [7] - - 


- 2J ( n » mn » + 0- - n i) 

i 


In(l-nj)) 


(44) 



We cannot just use a diagrammatic expression for 
from FT-MBPT as a functional in FT-RDMFT because 
in the framework of FT-MBPT, the noninteracting and 
interacting systems have different eq-lRDMs. 

To overcome this conceptual problem, we introduce an 
effective, generally nonlocal, one-particle potential which 
counteracts the effect of the interaction on the eq-lRDM. 
This idea is similar to the approach of Gocrling and Levy 
[45] in the context of zero-temperature DFT. By employ- 
ing a coupling constant integration we are then able to 
relate the interacting system to a noninteracting system 
with the same eq-lRDM. 

The grand potential functional of the fully interacting 
system reads 



n[ 7 ] = 

Tr{%] ( 



T 



+ V ext -liN + W-l/p\nD[ 7 ])}. (45) 



D(X) now represents the eq-SDO corresponding to the 
Hamiltonian f + XW + V(X). Using the fact that V(X) is 
a one-particle operator, the grand potential functional 
evaluated at the eq-lRDM -f eq (x,x') can now be ex- 
pressed as 



0[ 7 , 



<-</J 



0(1) 

n(o) 



1 dfi(A) 
~dX~ 



(49) 
(50) 



As 0(A) describes a system in thermal equilibrium it 
can be written with the help of the partition function 
Z(X) = t r { e -P(T+\w+v{\))} ag q( A ) = -i/^in^A). 

Furthermore, the equilibrium property implies that a 
variation of D(X) alone leaves 0(A) invariant. Taking 
the derivative in Eq. (50) w.r.t. the coupling constant 
therefore yields only contributions from XW and t^(A). 



dO(A) 
dX 



= Tr <^ D(X) W 




(51) 



V(X) is a one-particle operator and, since by construc- 
tion the eq-lRDMs of the interacting and noninteracting 
systems are equal for all A, Eq. (51) becomes 



rfO(A) 
dX 



= tr{D(X)W} 



,dv{X){x',x) , 

dxdx j(x, x ). (52) 

dX 



Using this result in Eq. (50) gives us our final expression 
for the grand potential functional evaluated at the eq- 
1RDM. 



0[7 eg ] = O fc [7 e9 ] + V ext [jeq} ~ pN[) eq \ 



l//?So[7e 



dATr 



|£>(A)#| (53) 
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Comparing Eq. (53) with Eq. (40) yields our final ex- 
pression for f2v^[7] as a coupling constant integral. 



Slwb] =Tr{£>{\)w} 



d\(W}> 



(54) 



In order to construct approximate functionals for the 
grand potential we need to evaluate J Q d\(W)\ which 
can be done by means of FT-MBPT. 
It is known (e.g. [46], p. 230) that (XW)\ can be related 
to the Green's function Q x of the system. 



1 



(XW)x = - / dxdx' lim - 5{x - x') 

2 J T'— >T+ \ Ul 

- (- S(x - x')^- + v{X){x' \x)^jg x {xT,x l T l ). (55) 

This expression, however, is not suitable for an imple- 
mentation into FT-RDMFT, because the noninteracting 
system, governed by H + V(X) exhibits a different eq- 
IRDM than the interacting one defined by Hq + XW + 
V"(A). Consequently the noninteracting Green's functions 
in a perturbativc expansion would not be easily deduce- 
able from the 1RDM which was given as an input to the 
functional. We therefore have to work with the nonin- 
teracting system defined by H + V(0) as our reference 
system in a perturbative expansion. Thus, Eq. (55) is 
rewritten as 



(XW)) 



dxdx' lim 



5{x — x ) 



- ( - 6(x -arO^- + v(0)(x', x))jg x (xr, x'r') 

+ i J dxdx'(v(0)(x', x) - v(X)(x', x))i(x, x'). (56) 

Now, the well known methods of FT-MBPT can be ap- 
plied to this unperturbed Hamiltonian. The perturbation 
consists of a two-particle interaction XW and a (possibly 
nonlocal) one-particle potential U(X). 



u{X) = v(X) - v(0) 



(57) 



The proof of Wicks theorem is still applicable for this 
kind of perturbation and the same Feynman rules apply. 
We show our notation conventions in Table I. 

As a short excourse, we would like to point out that if 
the Hamiltonian is temperature-independent and if the 
system is uniform then Eq. (56) can be written entirely 
in terms of Feynman diagrams. 



(XW)l nif = 




(58) 





Xw(x, x') 


vwvww 


u(X)(x, x') 


> 


g°(x,r,x',r') 




g x (x,T,x',T') 



TABLE I: Notation conventions for Feynman diagrams in FT- 
RDMFT 



Going back to general systems, the irreducible self en- 
ergy S* for the first-order contribution becomes 



E* = 



(59) 



Combining Eqs. (54), (56) and (59), we arrive at the 
first-order contribution to the interaction-induced grand 
potential functional in FT-RDMFT. 



^h[j) — 2 / dxdx'w(x, x') r y(x, x) r y(x' , x') 
Q x [y] = — ^ J dxdx'w(x, x') r y(x, x') , y(x', x). 



(60) 
(61) 
(62) 



where S* denotes the irreducible self-energy. 



We see that the functional form of the first-order con- 
tributions are equivalent to the Hartree- and Fock- 
functionals in zero-temperature RDMFT. Accordingly, 
we denote the remainder as correlation and get our fi- 
nal expression for the grand potential functional in FT- 
RDMFT. 



fi[ 7 ] = OfcH + Vexth] - A* AT [7] - l//?SoH 

+ fiff[7]+n x [ 7 ] + ft c [7] (63) 

Approximations for the correlation functional can now 
be derived by expanding the Green's function to higher 
orders and then solving Eqs. (54) and (56). 

We show in the following that the denotation as 
Hartree- and Fock-functionals is appropriate because the 
minimization of the grand potential using the function- 
als (60)-(62) is equivalent to a solution of the finite- 
temperature Hartree-Fock equations. 
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D. Finite-temperature Hartree-Fock (FT-HF) 

We consider the functional from Eq. (63) without the 
correlation part: 



n HF [ 1 } = n k [ 1 ] + v ext [ 1 ] 



fiNfr] - l//3S [l] 
+ fiir [7] + fix bl- 



The derivative of the Kohn-Sham entropy So [7] w.r.t. the 
occupation numbers is divergent for m — > {0, 1} whereas 
all other contributions are finite. Therefore, at the min- 
imum of the functional ^#^[7] there will be no pinned 
states. Furthermore, Eq. (64) constitutes an explicit 
functional of the 1RDM. We conclude that the functional 
derivative w.r.t. the 1RDM exists and that at the mini- 
mum the functional fulfills the following Eulcr-Lagrangc 
equation. 



gfigF[7] 
6j(x', x) 



= 



(65) 



We will now apply this condition on fi^f^M from Eq.(64) 
and project the result on the i-th natural orbital of the 
1RDM. This then leads to the following FT-HF equa- 
tions. 



o=/<**W)^ [7] 



V 2 



&y(x', x) 
4>i(x) + / dx'v ext (x,x')4>i(x') 



(66) 



dx'w(x, x') r y(x, x')4>i (x ! ) 



dx'w(x,x') r )(x' ,x')4>i(x) — Si4>i, (67) 



where in the last term we used Eq.(18). 

The physical situation whose equilibrium is found by 
a minimization of the grand potential implies the cou- 
pling to a particle bath. There are, however, important 
physical problems in which the particle number is intrin- 
sically conserved and therefore the equilibrium is defined 
by the state which minimizes the free energy instead. 
An important representative of these problems is, e.g., a 
molecule in solution where the solvent is described by a 
bath of harmonic oscillators at finite temperature. We 
will therefore conclude our investigation of FT-RDMFT 
by considering canonical ensembles in the following. 



E. Canonical ensembles 

The appropriate Hilbert space for the description of 
canonical ensembles of particle number N is given by the 
iV-particle subspace H N of the Fock space TL of Eq. (2) 



u 



N 



Sh 



(68) 



and the SDOs are weighted sums of projection operators 
on H N : 



W aN > 0, 2J WaN = 1. (69) 



The variational principle which governs this situation in- 
volves the free energy F = E — 1//3S rather than the 
grand potential and the eq-SDO is given by 



-/3H 



d: 



tr{e-P H }' 



(70) 



where H is now the A-particlc Hamiltonian of the sys- 
tem. 

Considering the framework of FT-RDMFT, the main 
difference to the grand canonical ensemble occurs when 
we investigate the canonical Kohn-Sham system. As in 
the grand canonical case the NOs of an eq-lRDM of a 
nonintcracting system will be given by the eigenstates of 
the one-particle Hamiltonian and the ONs will lie in the 
interior of the set of ensemble N-represcntablc lRDMs. 
However, because there is no simple analytic relation be- 
tween the eigenvalues of the one-particle Hamiltonian 
and the ONs as in Eq. (17), we do not know if every 
1RDM with < rii < 1 is a canonical eq-lRDM. If we 
know that a given 1RDM corresponds to a canonical equi- 
librium the corresponding potential has to be found by 
iterative methods similar to [35]. 

In our investigation of the Hohenberg-Kohn theorems 
for grand canonical ensembles we already mentioned that 
the one-to-one mapping between the eq-SDO and the eq- 
1RDM stays valid also in the case of canonical ensembles 
with the only difference that the external potential is now 
determined only up to an additional constant. Therefore, 
one can define a canonical universal functional 5" c [7] on 
the whole domain of ensemble-N-representable lRDMs 
as 



ben N - 



tT{D(f + W + l/pinD)}. (71) 



The equilibrium of the system is then found by a min- 
imization of the free energy functional F[y] = S c [l] + 

v ext [i\. 

So far, the formulation of FT-RDMFT for canonical 
ensembles follows in the same steps as for grand canon- 
ical ones and is equally well founded. Only the Kohn- 
Sham potential corresponding to a given 1RDM cannot 
be constructed as easily as in the grand canonical case. 
The reason for this lies in the fact that the exponent of 
the eq-SDO in a canonical ensemble contains only TV- 
particle contributions. Unfortunately, the finite temper- 
ature version of Wick's theorem [47] cxplicitely relies on 
the interplay of states of different particle numbers and 
is therefore invalid for canonical ensembles. Accordingly 
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the perturbative methodology as derived by us in the pre- 
vious section can generally not be translated to canoni- 
cal ensembles. However, considering the special case of 
the system being in the thermodynamic limit where the 
thermodynamic variables, and therefore also the corre- 
sponding functionals, of grand canonical and canonical 
ensembles coincide we can use a perturbative functional 
for the grand potential 51 [7] to calculate the free energy 

F[7] = n[ 7 ]+ M JV[ 7 ]. (72) 

IV. SUMMARY AND CONCLUSIONS 

In this work, we have shown that the equilibrium prop- 
erties of a grand canonical ensemble with local or non- 
local external potential are determined unambiguously 
by the eq-lRDM. We have used this fact to establish a 
functional theory for the grand potential in terms of the 
1RDM where the domain is defined easily in terms of sim- 
ple constraints on the occupation numbers of the 1RDM. 
The functional can then be divided in a universal part, 
independent of the external potential, and a part solely 
describing the external potential contribution. The eq- 
1RDM is then accessible by a minimization of the grand 
potential functional. 

Furthermore, we have shown that in contrast to the 
zero-temperature case, at finite temperature there ex- 
ists a nonintcracting Kohn-Sham system, reproducing 
the eq-lRDM of the interacting system under consid- 
eration. The existence of this effective nonintcracting 
system opens up several opportunities for further the- 
oretical development. Firstly, it has allowed us to de- 
velop a methodology for the construction of function- 
als in FT-RDMFT, using methods from FT-MBPT. Sec- 
ondly, one can incorporate the Kohn-Sham system in a 
self-consistent minimization scheme for FT-RDMFT, a 
project of active research by the authors. And finally, 
the Kohn-Sham system might grant access to single par- 
ticle properties which could be used for interpretation of 
the many-body properties of the system. 

As the resulting first-order functional is capable of de- 
scribing the kinetic energy and the exchange energy ex- 
actly, one only has to approximate the contributions from 
the interaction and the entropy. In comparison with the 
success of RDMFT for zero temperature, this sparks the 
hope that problems, which were not successfully treated 
by FT-DFT, could be accessed more easily in the theo- 
retical framework of FT-RDMFT. Furthermore, we have 
shown that the minimization of the first-order functional 
is equivalent to the solution of the finite-temperature 
Hartree-Fock equations. 

Finally we have considered the description of a quan- 
tum system in canonical equilibrium and have discussed 
how the results from the grand canonical ensemble differ. 
We have found that the description of canonical ensem- 
bles via FT-RDMFT is equally well founded. Its appli- 
cation however is hampered by the loss of the direct an- 
alytic relation between the Kohn-Sham Hamiltonian and 



the corresponding 1RDM and the invalidity of Wick's 
theorem for canonical ensembles. Nonetheless, under the 
assumption of being in the thermodynamic limit, func- 
tionals of FT-RDMFT for grand canonical ensembles can 
be used to calculate the thermodynamic variables of a 
corresponding canonical ensemble. 

After the theoretical foundation of FT-RDMFT, as 
presented in this article, the main task for the future is 
the development of correlation functionals for the grand 
potential and free energy in FT-RDMFT and the ap- 
plication to real systems. Some aspects of this will be 
discussed in parts II [27] and III [28] of this series. 

Appendix A: Equilibrium occupation numbers in 
general systems 

As we have pointed out in Section III A, the eq-lRDM 
of a noninteracting system will have occupation numbers 
strictly between and 1. We will now show that this 
is also true for the occupation numbers of eq-lRDMs of 
arbitrary systems, including interacting ones. 

We start from the spectral representation of the eq- 
1RDM. 

1 {x,x') = Y,ni4>* i {x')<l>i{x) (Al) 

i 

The occupation number operator hi is now defined as 

h, = c+c, (A2) 

where cf (cj) creates (annihilates) the natural orbital 0j. 

An arbitrary occupation number of the eq-lRDM in 
grand canonical equilibrium can then be written as 

m =tr{Dhi} =J2 w e(^e\h t \^e)- (A3) 

e 

The {^ e } are eigenfunctions of the Hamiltonian of the 
system and form a basis of the underlying Hilbert space. 
Another basis is formed by the Slater determinants {$ Q } 
which are constructed out of the natural orbitals {<fii} of 
the eq-lRDM of the system. The transformation between 
these bases is governed by the expansion coefficients c ea 
via 

* e = ^c eQ $ Q . (A4) 

Due to completeness and normalization of the {^e} and 
{$ Q }, the coefficients fulfill 

5> eQ | 2 = 5> eQ | 2 = l. (A5) 

e a 

Expanding the {^ e } in Eq. (A3) in terms of the {§ a } 
then leads to 

m =^w e ^c* ea c ef i{^ a \h l \^p) (A6) 

e a/3 
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Since the Slater determinants {$ Q } are constructed to 
be cigcnfunctions of the occupation number operator n,, 
this reduces to 



EE 



W e C e 



l*a) 



(A7) 



Using Eq. (A5) and the properties of the thermal 
weights, w e > and J2 e w e = 1; we see t na * 



fa >o 



(A8) 
(A9) 



The factors gi a are equal to 1 if the natural orbital 4n 
appears in the Slater determinant $ Q and otherwise. 
The summation over a corresponds to a summation over 
a basis of the Hilbert space, which is in the situation of a 
grand canonical ensemble the Fock space. Therefore, for 
a fixed i, there will be at least one a, such that gi a = 1 
and at least one a for which gi a = 0. Applying this 
result, together with Eqs. (A8) and (A9), to Eq. (A7) 
yields the desired inequality relations. 



< Hi < 1 



(A10) 



Appendix B: Zero-temperature mapping between 
potentials and wavefunctions 

To investigate the nature of the mapping between the 
nonlocal external potential and the corresponding gs- 
wavcfunction we will consider the relationship between 
the potential and the gs-lRDM and then use Gilbert's 
theorem [26] to translate the results. 

Assume an arbitrary Hamiltonian H with gs-lRDM 



H = T + V + W 



(Bl) 
(B2) 



Due to Gilbert's theorem, the wavefunction can be writ- 
ten as a functional of the 1RDM and therefore an ener- 
gyfunctional E[j] can be defined. 



s[ 7 ] = {ni]\H\m\) 



(B3) 



From the variational principle, this functional is minimal 
for the gs-lRDM and therefore the following relations 
hold 



dm 



a, > /i , m = 

fj, , < m < 1 

bi < fx , rii = 1 



(B4) 



where \i can be identified as the chemical potential of 
the system. We will now be able to show that the map- 
ping between potential and gs-lRDM is one-to-one if and 
only if there are no pinned occupation numbers, i.e. oc- 
cupation numbers equal to or 1. This will be done in 
two steps. At first we show that in the case of unpinned 
occupation numbers the external potential is uniquely 
determined up to a constant. Secondly, we will consider 
gs-lRDMs with pinned occupation numbers and show 
that one can explicitly construct infinitely many poten- 
tials which leave the gs-lRDM invariant. 

• unpinned states: 
The absence of pinned states allows us to use the follow- 
ing Euler-Lagrange equation 



6E[>y] 



Sj(x, x 1 ) 



5N[<y} 



5j(x, x') 



fid(x, x 1 ) 



(B5) 



7gs 



The addition of an arbitrary potential contribution U[j] 
to the energy functional E u [y] = E[y] + U[j) will then 
yield 



SE u [>y] 



5j(x, x') 



fxS(x, x 1 ) + u{x, x') 



(B6) 



As every gs-lRDM has to fulfill Eq. (B5) we can deduce 
that the only possible choice of u(x, x') which leaves the 
gs-lRDM invariant is 



u(x, x') = cS(x, x'). 



(B7) 



• pinned states: 
As in the case of pinned occupation numbers the min- 
imum of E[y) will be at the boundary of the domain, 
we cannot use Eq. (B5). It would be possible to ad- 
just the Euler-Lagrange equation by an incorporation of 
Kuhn- Tucker multipliers [41] but there is a simpler way, 
as described in the following. 

The fact that for pinned states, the derivatives are al- 
lowed to be different from /i allows us to construct one- 
particle potential which leave the gs-lRDM invariant. 
This potential shall be governed by the generally non- 
local kernel u(x, x'). By choosing it to be diagonal in the 
basis of natural orbitals of the gs-lRDM we ensure that 
the orbitals will not change upon addition of the poten- 
tial. For simplicity, we choose only one component to be 
non- vanishing. 



u(x,x') = u(j)* a (x')4> a (x) 
The new energyfunctional is then given by 

-E-aM = E[j] + / dxdx'u(x,x / )j(x' ,x) 



(B8) 



(B9) 



and the derivative with respect to the occupation num- 
bers becomes 



dE a [^} 



drii 



CH + u5ia , Ui = 

jJL + uSia , < Hi < 1 

bi + uS ia , Hi = 1 



(BIO) 
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These considerations can now be employed to show 
the ambiguity of the external potential in RDMFT for 
groundstates with pinned occupation numbers. For sim- 
plicity we will now assume that there is exactly one 
pinned occupation number with np = 0. We then con- 
struct a potential of the form of Eq. (B8) with a = 0. 
From Eq. (BIO) we see that every choice of u > /i — ap 
leads to a situation where the /3-orbital exhibits a deriva- 



tive bigger than /i, leaving the gs-lRDM invariant. For 
one pinned ocuppation number np = 1 we can choose 
u < [i — bp which lets the derivative of the /3-orbital al- 
ways surpass n which again leads to the same gs-lRDM. 
When considering gs-lRDM with several pinned states 
these arguments are readily translated which proofs the 
ambiguity of the one-particle potential for gs-lRDM with 
pinned occupation numbers. 
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